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ABSTRACT 

We suggest that Jovian planets will survive the late stages of stellar evolution, and 
that white dwarfs will retain planetary systems in wide orbits (> 5AU). Utilising 
evolutionary models for Jovian planets, we show that infra-red imaging with 8m class 
telescopes of suitable nearby white dwarfs should allow us to resolve and detect com- 
panions > 3 Mjup . Detection of massive planetary companions to nearby white dwarfs 
would prove that such objects can survive the final stages of stellar evolution, place 
constraints on the frequency of main sequence stars with planetary systems dynami- 
cally similar to our own and allow direct spectroscopic investigation of their composi- 
tion and structure. 
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1 INTRODUCTION 

Over 70 extra-solar planets have now been detected since 
the discovery of a companion to the solar-type star 51 Peg 
in 1995 (Mayor & Queloz 1995). All of these planets have 
been discovered via the radial velocity technique, in which 
the planet's presence is inferred by the central star's mo- 
tion around the system's barycentre. Instrumental and in- 
trinsic noise limit the sensitivity of these studies to Saturn- 
mass companions in short-period orbits, with more massive 
companions detectable in more distant orbits. The dura- 
tion of current programmes limits this to ~3.5 AU (Fischer 
et al. 2002), although some systems do display interesting 
trends indicative of more distant companions. Therefore, 
we have little information on systems with massive com- 
panions at large radii, such as our own Solar System. The 
only constraint on Jovian-like systems comes from micro- 
lensing statistics, which suggest less than 1/3 of M stars 
have Jupiter-like planets orbiting at 3 — 7 AU ( Albrow et al. 
2001). 

The planets discovered by the radial velocity technique 
are not open to further direct study, due to their close 
proximity to the much brighter parent star. The only ex- 
ception to this is the transiting companion to HD209458 
(forb ~ 3.5 days). Transit photometry shows this planet is 
a gas-giant, and Charbonneau et al. (2002) have recently 
detected sodium in its atmosphere. Still, the planet itself 
cannot be directly imaged. 

Several groups have conducted imaging surveys of 
nearby main sequence stars to search for low mass com- 
panions, including the use of adaptive optics, coronographs 
and space-based observations with HST (e.g. Turner et al. 



2001, Neuhauser et al. 2001, Oppenheimer et al. 2001, Kuch- 
ner & Brown 2000, Schroeder et al. 2000). However, the 
extreme contrast (> 20 magnitudes) and small separations 
(5AU= 1" at 5pc) between main sequence stars and Jovian 
planets makes sensitive surveys very difficult. To date, no 
planetary mass companions to nearby stars have been im- 
aged, although several brown dwarf companions have been 
detected via direct imaging (e.g. Nakajima et al. 1995). 

The end state of main sequence stars with M<8M0, 
white dwarfs, are typically lO'^-lO* times less luminous than 
their progenitors. Thus, there is potentially a strong gain in 
the brightness contrast between a planet and a white dwarf 
when compared to a main sequence star, assuming that 
planets can survive the late stages of stellar evolution. The 
gain in contrast is strongest in the mid-infrared, where the 
planet's thermal emission peaks well into the white dwarf's 
Rayleigh- Jeans tail. 

Indeed, Ignace (2001) has suggested that excess infra- 
red emission could be detectable from a hot Jupiter orbiting 
a 10,000K white dwarf at a distance of ~ 10"^ white dwarf 
radii and with an orbital period of ~ 10 days. Chu et al. 
(2001) have also suggested that, near a hot UV-bright white 
dwarf, the atmosphere of a Jovian planet would be photoion- 
ized and emit variable hydrogen recombination lines, which 
may be detected by high-dispersion spectroscopic observa- 
tions. However, both these methods rely on the planet being 
in a close (0.01-2AU) proximity to the white dwarf, where 
it would be difficult to resolve. 

In this letter, we discuss the potential for imaging plan- 
etary companions in wide (>5AU) orbits around nearby 
white dwarfs. In section 2 we investigate the probability of 
a planetary system surviving the late stages of stellar evolu- 
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tion. We discuss the dotectability of any surviving planets in 
section 3, suggest suitable target white dwarfs in section 4, 
and discuss the expected frequency of wide planetary com- 
panions to white dwarfs in section 5. 



2 PLANETARY SYSTEMS IN THE 
POST-MAIN SEQUENCE PHASE 

Any planetary companion to a white dwarf must have sur- 
vived the red giant branch (RGB) and asymptotic giant 
branch (AGB) phases of stellar evolution. During these 
stages, the white dwarf progenitor swells up to a few hun- 
dred solar radii in size and undergoes significant mass loss. 
What happens to a planet during this time depends on the 
initial orbital separation, the stellar mass loss rate, the to- 
tal mass lost, tidal forces, and interaction with the ejected 
material. 

Any planet in an initial orbit within the final extent 
of the red giant's envelope will be engulfed and migrate in- 
wards. Livio & Soker (1984) showed that these planets will 
either completely evaporate or accrete mass and become a 
close companion to the eventual white dwarf. 

Planets in wider orbits that avoid direct contact with 
the expanding red giant envelope (such as the Jovian plan- 
ets of our own solar system) will have a greater chance of 
survival, migrating outward as mass is lost from the central 
star. However, stars more massive than the Sun lose more 
than half of their original mass in evolving to the white dwarf 
stage (Weidemann 1987, Wood 1992). If over half the mass 
of the central star is lost instantaneously then any accompa- 
nying planets will likely escape, since their orbital velocity 
will suddenly exceed the escape velocity of the system. Thus, 
we might naively expect no white dwarfs to retain planetary 
systems. 

Of course, in reality, main sequence stars of M<8.0 M0 
do not evolve into white dwarfs instantaneously. Mass loss 
on the RGB and AGB, and final ejection of the planetary 
nebula is a relatively slow process. The Sun, for example, 
will lose around 0.2 M© on the RGB in 10® years, and 
up to another 0.2 M0 on the AGB on a timescale of up to 
~ 10^ years. The final thermal pulse, in which the cool red 
giant develops a super-wind and finally ejects the remainder 
of its envelope, exposing the core, is poorly understood but 
probably occurs on a timescale of a few x 10* years (Soker 
1994). 

The dynamical timescale for a planet to react to the 
change in the mass of the central star, and thus the gravi- 
tational force between the star and planet, is given by; 
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where a is the semi-major axis of the orbit and Vorb is the 
orbital velocity, and; 
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where Mi is the mass of the central star and M2 is the mass 
of the planet. Assuming Mi >> M2; 
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This is much less than the mass loss timescales men- 
tioned above. Even if a 2 M0 star lost its entire envelope 
(1.4 Mq) during a short FN phase of 10^ years (which it 
doesn't, only the last few tenths of a solar mass), the orbit 
would not become unbound. Thus, to a first approximation, 
the orbits of planets which do not interact directly with the 
red giant will simply expand adiabatically by a maximum 
factor Mms/Mwd (Jeans 1924, Zuckerman & BeckUn 1987). 
For example, for a main sequence star ^ 2 Mq and a white 
dwarf ~ 0.6 Mq, the orbit would expand by a maximum fac- 
tor ~ 3. A planet in a Jupiter-like orbit would be relocated 
to ~ 15 AU, and a planet in a Neptune-like orbit to ~ 90 AU. 

For planets in orbits <10AU, tidal forces between the 
planet and red giant are important in slowing orbital mi- 
gration by transfering angular momentum from the orbit to 
stellar rotation. For example, Soker (1994) shows that, ne- 
glecting tidal forces, as the Sun evolves to a white dwarf 
Jupiter's orbital radius will expand from its present value 
of a = 1118-Rq to a = 1860i?Q. However, tidal forces on 
the RGB will counteract this expansion, reducing this final 
orbital radius by ~ 3%, and by a further ~ 14% on the 
upper AGB. For planets within 5 AU, Livio & Soker (1983) 
showed that tidal interactions will actually cause them to 
migrate tnwards. Soker (1996) suggests that tidal interac- 
tions between evolving red giants and low mass companions 
could be responsible for producing the large fraction of el- 
liptical PN observed, indicating low mass companions may 
be common amoung PN progenitors. 

The orbits of distant planets will also be subjected 
to drag resulting from their interaction with the material 
ejected by the central star. Duncan & Lissauer (1998) de- 
fine the total amount of material impacting a planet during 
post-main sequence evolution as; 
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where Vsw, Vorb and Vesc are the velocities of the stellar wind, 
the planet and the escape velocity of the planet respectively, 
and Rp is the radius of the planet's magnetosphere. Evalu- 
ation of this equation, however, shows that gas giants will 
collide with significantly less than 1% of their own mass over 
the entire red giant phase. The orbits will decay inwards only 
slightly as a result. 

Duncan & Lissauer (1998) have made detailed simu- 
lations of the effects of post-main sequence mass loss on 
the stability of the solar system and for planetary systems 
around more massive stars. They show that the orbits of 
the giant planets will probably be stable for tens of billions 
of years subsequent to the Sun's death. However, planetary 
systems dynamically similar to our own around somewhat 
more massive stars may eventually be de-stabilised. Duncan 
& Lissauer's calculations show that giant planets around, for 
instance, an initially 4Mq star, will be liberated by large- 
scale chaos in less than a Hubble time. 

Finally, Soker (1999) points out that the newborn plan- 
etary nebula central star is hot (Teg > 100, OOOK) and a 
strong source of ionizing soft x-ray and UV radiation. The 
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interaction of this radiation with a planet may lead to ab- 
lation of the planetary atmosphere, but this will only be 
effective for planets having an escape velocity Vesc ^ Cs, 
where Cs = 15 kms~^ is the sound speed of the ionised plan- 
etary material. For a Jupiter-mass planet Vesc = 61kms~^ 
and thus there will be negligible ablation. 

It seems likely, then, that giant planets in initially dis- 
tant (a>5AU) orbits around main sequence stars 1 — 8 M0 
will survive the late stages of stellar evolution, and remain in 
orbit around the remnant white dwarfs at increased orbital 
radii for at least several billion years. 



3 DETECTING PLANETS AROUND NEARBY 
WHITE DWARFS 

Our ability to detect an extra-solar planet around a nearby 
white dwarf will depend on its intrinsic luminosity, which 
in turn depends on its age, its distance from us, and its 
separation from the parent star. 

3.1 Luminosities of extra-solar planets 

Burrows et al. (1997) have made nongray calculations of the 
expected spectra, colours and evolution of solar-metallicity 
massive planets. We have used these models to predict the 
infra-red magnitudes of planets around nearby white dwarfs. 
The Burrows et al. models describe the evolution of a planet 
in isolation, and do not include the effects of thermal insola- 
tion by, or reflected light from, a parent star. In reality, gas- 
giant planets will have been kept warm throughout their 
lives by the white dwarf progenitors. They will also have 
accreted matter during the RGB/AGB/PN phases. Both ef- 
fects would act to increase their temperatures and luminosi- 
ties, although in the white dwarf phase reflected light makes 
a negligible contribution to the planet's brightness. Our pre- 
dicted luminosities are, therefore, lower limits. 

3.2 The ages of white dwarfs 

The white dwarf cooling age and progenitor main sequence 
lifetime are the dominant timescales determining the bright- 
ness of any planetary companions, as they are significantly 
longer than the RGB or pre-main sequence phases. 

White dwarf cooling ages are dependent on a number 
of parameters including: effective temperature, mass (usu- 
ally derived from the surface gravity), composition of the 
core and the thickness of the surface H- (or He-) dominated 
atmosphere. Temperature, surface gravity and atmospheric 
composition are usually determined through fitting model 
atmospheres to optical and UV spectra (e.g. Bergeron, Saffer 
& Liebert 1992, Marsh et al. 1997, Finley, Koester & Basri 
1997). The mass is then inferred from the surface gravity 
measurement using evolutionary models (e.g. Wood 1992, 
1995). These same models also give the cooling age. The 
cooling ages and masses of many nearby white dwarfs have 
been calculated in this manner by Bergeron, Leggett & Ruiz 
(2001, BLR). For nearby white dwarfs not included in the 
BLR sample (generally, those hotter than ~ 10, 000 K), we 
have estimated cooling ages from published temperature and 
meiss measurements. 

The relationship between the mass of a white dwarf 



and its main sequence progenitor, the initial-final mass rela- 
tion, has been investigated by a variety of authors. We have 
adopted model A from Wood (1992): 
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Wood (1992) also gives a simple formula for estimating main 
sequence lifetimes: 
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The white dwarf cooling age and the main sequence lifetime 
can then be combined to give the system age of local white 
dwarfs. We use this age, combined with the distance to the 
white dwarf, to estimate the magnitude of planets. 



4 SUITABLE TARGETS AMONG NEARBY 
WHITE DWARFS 

There are 118 white dwarfs currently catalogued within 20pc 
of the Sun (Holberg et al. 2001). All of these systems are 
close enough to allow us to resolve planetary companions in 
^ 5-100 AU orbits from the ground. For example, at 10 pc 
an orbit of 10 AU is equivalent to a maximum separation of 
1" on the sky. Planets in wider orbits will be even easier to 
resolve. Adaptive optics, coronographs and techniques such 
as nulling interferometry will increase the chances of resolv- 
ing closer companions. 

However, only those white dwarfs young enough for 
Jovian companions to still be detectable will make suit- 
able targets for extra-solar planet searches. We are most 
likely to be able to detect an extra-solar planet around 
a relatively warm and young (<a fewxlO^yrs), massive 
(>O.6M0) white dwarf descended from a relatively mas- 
sive progenitor (>2M0), since its main sequence lifetime 
is short (<1.8Gyr). For Mwd<O.6M0, the main sequence 
lifetime is very sensitive to errors in the Mwd estimate (see 
Figure |l|). Therefore, we would be unwise to rule out white 
dwarfs with mass estimates below 0.6 M0, as they may be 
descended from more massive stars than we expect. 

Table |l] lists the top 10 candidates among nearby white 
dwarfs for a search for planetary companions, considering 
the main sequence lifetime Ims, the white dwarf cooling age 
tcooL and the distance. Also included in the table are the 
estimated J band magnitude of 3, 5, & 10 Mjup planets in 
each system. 

In Figure ^ we plot the magnitude of a 5M,iup 
planet around these white dwarfs in several near-IR bands 
(J,H,K,L' & M). Also plotted are the nominal sensitivity 
limits for current near-infrared imaging instruments (from 
on-line exposure time estimates); Gemini-North -I- NIRI 
(solid lines) and VLT + ISAAC (dashed lines). We can see 
that although the planets are brighter in the M band {5fim), 
the instruments are only sensitive enough in the J and H 
bands (1.1 & 1.6/im). 
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Table 1. Top 10 candidate white dwarfs for detecting Jovian planets. White dwarf parameters are taken from; Bergeron, Saffer & 
Liebert (1992; BSL92); Wickramasinghe & Ferrario (2000; WFOO); Bergeron et al. (1995; B95); Bergeron, Leggett & Ruiz (2001; BLROl); 
Bragaglia, Renzini & Bergeron (1995; BRB95) and McMahan (1989; M89). 
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Figure 1. The relationship between white dwarf mass and main- 
sequence lifetime (from equations ^ and The steepness of the 
relationship in the range 0.55<Mwd<0.6 shows why we should 
not rule out white dwarfs with mass estimates in this range, de- 
spite the indications their progenitors had a long main-sequence 
Ufetime. 




Figure 2. Magnitudes of a 5Mjup planet around nearby white 
dwarfs (Table hi) in the near-IR J,H,K,L' and M bands. Also plot- 
ted are the nominal limiting sensitivities (Ihr, b-a) for Gemini- 
N-I-NIRI (solid lines) and VLT+ISAAC (dashed lines). 
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5 THE EXPECTED FREQUENCY OF 

PLANETS IN WIDE ORBITS AROUND 
WHITE DWARFS 

Assuming that massive planets born in initially wide (> 
5AU) orbits remain bound to white dwarfs, then the ex- 
pected frequency of massive planets around nearby white 
dwarfs is simply a function of the frequency of solar-like and 
more massive main sequence stars with Jovian-like plane- 
tary systems. Unfortunately, we currently have little obser- 
vational information on these systems, although it is possi- 
ble to extrapolate from existing data to make some estimate 
of what we might expect. Lineweaver & Grether (2002, LG) 
have analysed the statistics of planets detected via radial ve- 
locity techniques, and estimated the number of Jovian-like 
systems we might expect to find. They have produced a fit to 
the completeness-correct period distribution of planets so far 
detected, in the mass range 0.84 Mjup < M sini <13M,jup; 
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= a log P + b 



(6) 



dlogP 

where a = 12 ± 3 and b — —9 ± 2. The initial orbital periods 
of interest range from ~ 10 yr, to avoid direct interaction 
with the expanding red giant envelope, to ~ 200 yr, where 
we might expect Jovian planet formation to end. Integrating 
equation]^ over this range, we estimate there to be ~55±16 
Jovian planets per ~1000 white dwarf systems. However, 
for a significant sample of nearby white dwarfs we can only 
detect planets > 5 Mjup with current technology. LG fit the 
observed planetary mass function as; 

diV 

— — — — — =alog(Afsini)-h6 (7) 
dlog(A/ smi) t,\ / \ / 

where a — — 24±4. To produce 55±16 planets over the mass 
range 0.84 Mjup < M sini < 13 Mjup we need to set b = 
— 58.4± 19.7. Therefore, from equation ^, we expect ~15±6 
planets > 5M jup per 1000 white dwarfs. Hence, we may 
expect ~l-2% of white dwarfs to possess one or more planets 
> 5 Mjup . Note that we treat Msini as equivalent to M, so 
this predicted number density is a lower limit. 

This analysis is, however, based on statistics for plan- 
ets in orbits < 3.5 AU. It is quite possible that laws derived 
for planets in this area of orbital parameter space around 
~1Mq stars do not apply to Jovian planets at large radii 
(>5 AU) around more massive stars (2 — 8M0). 
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6 DISCUSSION 

We have suggested that Jovian planets wiU survive the 
late stages of stellar evolution, and that some nearby white 
dwarfs possess planetary systems in wide orbits. Utilising 
evolutionary models for Jovian planets, we have shown that 
infra-red imaging of suitable nearby white dwarfs should al- 
low us to resolve and detect companions M>3 Mjup with 8m 
class telescopes. The best candidates for an observational 
search are relatively hot (young), massive (> 0.6 M©) white 
dwarfs, as they have the shortest overall system age. How- 
ever, given that we have neglected potential heating mecha- 
nisms such as thermal insolation and accretion from the red 
giant's wind, Jovian planets around nearby white dwarfs 
are likely to be brighter than we have suggested in Figure ^ 
and may be detectable around ~ 0.55 M© white dwarfs de- 
scended from progenitors as low in mass as ~ 1.2 Mq. 

Our "top ten" list of targets among nearby white dwarfs 
(Table 1) includes degenerates descended from main se- 
quence progenitors from 2.6 M© — 7.4 Mq, corresponding 
to spectral types early A and B. The existence of a dust 
disk around Vega (AOV) strongly suggests that planetary 
systems may exist around such early-type stars, although 
they are not the subject of current radial velocity searches. 

Detection of massive planetary companions to nearby 
white dwarfs would prove that such objects can survive the 
final stages of stellar evolution, place constraints on the fre- 
quency of main sequence stars with planetary systems dy- 
namically similar to our own, and allow direct spectroscopic 
investigation of their composition and structure. 

An unsuccessful search would suggest one of, or a combi- 
nation of, a number of scenarios, including: (a) that planets 
> 3 Mjup are fainter than expected from the evolutionary 
models we have used; (b) that, at least for a sub-sample 
of targets, any accompanying planets have too small a pro- 
jected separation from the white dwarf to be resolved; (c) 
that massive Jupiters in wide orbits around 1-8 M© main 
sequence stars are not common; (d) that these planets do 
not remain bound to white dwarfs. This last scenario would 
raise the possibilty that the Galaxy contains a population of 
isolated Jovian-mass planets ejected by white dwarfs during 
the final stages of their evolution. 
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